Turbulence Fluctuations and New Universal Realizability Conditions in Modelling 
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General turbulent mean statistics are shown to be charac- 
terized by a variational principle. The variational functionals, 
or "effective actions" , have experimental consequences for tur- 
bulence fluctuations and are subject to realizability conditions 
of positivity and convexity. An efficient Rayleigh-Ritz algo- 
rithm is available to calculate approximate effective actions 
within PDF closures. Examples are given for Navier-Stokes 
and for a 3-mode system of Lorenz. The new realizability 
conditions succeed at detecting a priori the poor predictions 
of PDF closures even when the classical 2nd-order moment 
realizability conditions are satisfied. 

PACS Numbers: 47.27.-i, 47.27.Sd, 47.11. +j 



It does not seem to be a well-recognized fact that gen- 
eral turbulence mean statistics — such as mean velocity or 
pressure profiles — are characterized by a variational prin- 
ciple. However, in nonequilibrium statistical mechanics it 
was pointed out long ago by Onsager |l],|| that mean his- 
tories satisfy a "principle of least action" . The so-called 
Onsager-Machlup action determines the probability of 
fluctuations away from the most probable state. Close 
to thermal equilibrium there is a standard fluctuation- 
dissipation relation, so that the action has the physi- 
cal interpretation of a "dissipation function" . Onsager's 
variational principle reduces then to a principle of least 
dissipation. 

Recently it has been pointed out by one of us 
that a similar effective action T[z] exists in turbulent flow 
for any random variable z(t). This action function (i) 
is non- negative , T[z] > 0, (ii) has the ensemble mean 
~z(t) as its unique minimum T[z] = 0, and (Hi) is con- 
vex, T[Xzx + (1 - \)z 2 ] > AT[2i] + (1 - A)T[z 2 ], < 
A < 1. These are realizability conditions which arise 
from positivity of the underlying statistical distributions. 
As a consequence, the mean value z(t) is character- 
ized by a "principle of least effective action". Just as 
Onsager's action, this functional is related to fluctua- 
tions. In particular, in statistically stationary turbu- 



lence, the timc-cxtcnsivc limit of the effective action, 
V[z] = Hindoo ~r[{z(i) = z : < t < T}}, the so-called 
effective potential, determines the probability of fluctu- 
ations in the empirical time-average z~t = ^ dt z(t) 
away from the (time-independent) ensemble-mean value 
z. More precisely, the probability for any value z of the 
time-average zt to occur is given by 



Prob ({z T « z}) ~ exp (-T • V[z\) 



(1) 



This agrees with the standard ergodic hypothesis, accord- 
ing to which, as T — > oo, the empirical time-average must 
converge to the ensemble-mean, z~t — ► z, with probabil- 
ity one in every flow realization. The Eq.([l]) refines that 
hypothesis, by giving an exponentially small estimate of 
the probability at a large (but finite) T to observe fluc- 
tuations away from the ensemble-mean. 

The realizability conditions on the effective action or 
effective potential hold even when there are no classical 
2nd-moment realizability conditions on the means them- 
selves. Thus, energy spectra or Reynolds stresses (2nd 
moments) must be positive, but mean velocity profiles 
(1st moments) or energy transfer (3rd moments) do not 
satisfy such simple realizability conditions. The new real- 
izability conditions thus have great potential significance 
in modelling if they can be efficiently calculated within 
turbulence closures. In j|M we have shown that there is 
a simple Rayleigh-Ritz algorithm within PDF closures — 
such as mapping closures ||U or generalized Langevin 
models |?]|| — by which the corresponding approximate 
values of the effective action or effective potential may 
be readily calculated. 

As a simple example, we consider first a 3-mode system 
of Lorenz M , in which the equations of motion are 



x>i Aix jXk v^xi F ft, 



(2) 



with i,j, k a cyclic permutation of 1,2,3, with A\ + 
A2 + As = imposed on interaction coefficients Ai 
for energy conservation, with vi positive damping coeffi- 
cients, and fi(t) white-noise random forces with covari- 
ance 2/tj. This dyamics has been used often as a first test 
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of closure ideas We consider a simple map- 

ping closure proposed by Bayly for the 3-mode system 
p3| , which models the realizations by a quadratic map 
Xi — /3iNi + PiN'jN'k of independent standard Gaussian 
variables Ni,N?, i — 1,2,3. The resulting closure equa- 
tions for the 2nd moments Mi = (xf), i = 1, 2,3 and the 
3rd moment T = {x\X2X 3 ) are 

Mi = 2A,T - 2viMi + 2m (3) 

for i = 1,2, 3 and 

T = AiM 2 M 3 + v4 2 MiM 3 + A 3 MiM 2 

-(^ + z/ 2 + j/ 3 )r. (4) 

These are just the quasinormal ( QN) equations for the 
3-mode system, obtained by neglecting the 4th-order cu- 
mulants |Q. It was already noted by Kraichnan |l(]] 
that, unlike for Navier-Stokes, the QN closure for the 3- 
mode system predicts all positive energies. In fact, for 
Ai = 2,A 2 = A 3 = -1,ki = 1,k 2 = k 3 = 0.001,z/i 
0.00 1 . v<x = U3 = 1 it gives steady-state values 

M[ QN) « 1.49875, M 2 (g7V) = Mp N) a 0.50025, 

y(QJV) ~ -0.49925. 1 j 

All of the 2nd moments are positive, as required by re- 
alizability. However, DNS of the 3-mode dynamics itself 
gives 

M (DNS) = 4 4g ± Q Q3 

m (dns) = M (QN) = o 49876 ± 0.00002 (6) 
rp(DNS) = _o.49776 ± 0.00002. 

While the QN predictions for M 2 ,M 3 and T are within 
j% of the DNS values, Mi is underpredicted by 66% 
in the QN approximation. As is well-known, satisfac- 
tion of realizability cannot guarantee that a prediction is 
correct. However, failure of realizability certainly implies 
that the predictions are in error. In Figs. 1-3 we graph the 
approximate effective potentials of the energy variables 
Ei — t^x\ and triple product T — xix 2 x 3 in the QN clo- 
sure as calculated by the Rayleigh-Ritz algorithm for the 
above PDF model. The numerical method is outlined be- 
low and described in detail in j|,[l5]]. It is apparent that 
Ve 2 and Vr satisfy realizability but that Ve 1 does not. 
Thus, one may infer a priori — without knowledge of the 
DNS results — that the QN prediction for the mean of E\ 
is not converged. In this case, the failure of realizabil- 
ity of the predicted Ve 1 succeeds at detecting the poor 
prediction for the mean value, even though the classical 
2nd-moment condition Ei > is satisfied. In the same 
plots in Figs. 2-3 we have graphed also the effective poten- 
tials Ve 2 and Vr obtained from DNS. They do not agree 
with the QN potentials as closely as do the corresponding 
means: the accurate prediction of fluctuations is a much 
more stringent demand on the closure. However, we note 



that the predictions of Bayly's PDF closure are at least 
qualitatively correct for Ve 2 and Vr and give correctly 
the order of magnitude of the averaging-time needed to 
eliminate fluctuations in those variables. Of course, the 
prediction of Ve 1 is not even qualitatively correct. 

The Rayleigh-Ritz algorithm used in obtaining the ap- 
proximate potentials from the PDF closure involves a 
fixed point problem very similar to (and, in fact, gen- 
eralizing) the fixed point condition determining the pre- 
dicted steady-state moments themselves. The system of 
equations that must be solved in general is 

|^(M,h)a +(^Q (M,h)-a = y (M,h)a, (7) 

V(M,h) = F (M,h)M, (8) 

a<o + a-M = 1. (9) 

The vector M = (Mi, ...,Mk) contains the moments of 
the closure, e.g., in the case above, k = 4 (and M4 = T). 
h is the vector of "perturbation fields" , one associated 
with each variable Zi for which the potential is to be de- 
termined. In our previous calculation h = (hE 1 ,h,E 2 i hr)- 
When h = the vector V(M, h) coincides with the dy- 
namical vector V(M) which appears in the closure equa- 
tion: M = V(M) (cf. Eqs.(§),(|) above). The perturba- 
tions for h are determined by the method discussed 
in 0]. The 0-component Vb(M, h) is associated with the 
zeroeth moment M = 1 and it may be written explicitly 
here as V (M,h) = \h El M x + \h E2 M 2 + h T M 4 . It is 
easy to check that, when h = 0, the stationary moments 
M* along with a*o = 1, a* = solve the system Eqs. ([?])- 
(^). Once the solutions a*o(h), ct*(h), M*(h) are known 
for h ^ 0, the effective potential Vz ( is constructed as a 
function of hi via VzJ/ij] = — a*(/ij)-V (M*(/ij)). To ob- 
tain the potential as a function of Zi, the expected value 
Z*(h) = z must be inverted to give hi as a function of 
Zi . For full details of the algorithm, see j||l5) . 

Our results point toward significant new directions in 
turbulence modelling. The new realizability conditions 
apply individually to all predicted means. We see above 
that they can successfully discriminate between poor pre- 
dictions for one set of variables and good predictions for 
another. Calculating each point on the graph of an ef- 
fective potential curve within a closure requires just the 
same amount of computation as that to calculate the 
predicted mean. It is therefore very easy to apply the 
above realizability conditions as a check to detect poor 
predictions in advance, without expensive testing by ex- 
periment or simulation. This gives a strong incentive to 
the development of PDF closures, such as those in Refs. 
p|- In conjunction with our variational method they 
can give some a priori information in turbulence mod- 
elling. This is a unique advantage, almost never obtained 
in other statistical closure methods. 
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It remains to be seen how well the new realizability 
conditions succeed in detecting poor predictions of clo- 
sures for Navier-Stokes turbulence. It is thus worthwhile 
to give one example demonstrating the Rayleigh-Ritz 
method for a statistically time-dependent Navier-Stokes 
flow. The simplest such situation is freely-decaying ho- 
mogeneous and isotropic turbulence with random initial 
data. We consider a model energy spectrum 



(p - 2)A m • (E(t) - K{t)) ■ KP-\t). 



(14) 



E(k,t) 



Ak m k<k L (t) 
ae 2 / 3 (i)fc" 5 / 3 k L (t) < k < k d (t) (10) 




k > k d {t) 



which has been adopted before in this problem Jl6| , p7t . 
As long as < m < 4 it is commonly believed that 
there is a permanence of the low-wavenumber spectrum. 
This motivates one to adopt the above self-preserving 
form, in which the shape of the spectrum is unchanged 
in time except through its dependence on the parameters 
(i) and kd(t). At high Reynolds number there is 
only one independent such parameter, since the relation 

3 

kh{t) = (% £ (t)) 3 '™ +5 is required by continuity and, when 

M*) « **(*)> k d (t) = (3^) 3/ V /4 W also holds @. 
The remaining time-dependence is determined by consid- 
ering the evolution of the mean energy E(t) — ^(v 2 (t)). 
For the above form of the spectrum it is not hard to show 
that the dissipation e(t) = %Y,ij(( d i v j + d j v i) 2 ) is 



given as 



e(t) = A m ■ EP(t) 



(11) 



3 \ 2m + 2 



, , v4™+i and p = t 2 ^. 

m + l 2 J * 2m+2 



with A" 1 = a 3/2 
Thus, employing the Navier-Stokes equation via its 
energy-balance, one obtains the closed equation 



E(t) = -A m • BP (t). 



(12) 



Its solution gives a prediction for the energy-decay law, 

as E*(t) ~ (t- toy^': see @. 

It is interesting to make a check on the various hy- 
potheses involved in these predictions by means of the 
effective action T[E] for the energy history E(t). As a 
simple PDF model for the above closure, one may adopt 
a Gaussian random velocity field with the assumed self- 
similar spectrum Eq.([To|). The Rayleigh-Ritz approxi- 
mation of the effective action within the Gaussian ansatz 
can be analytically evaluated with the result: 



y\(Gauss) rgj 

dt 



2(p - 2)A„ 

oo lE{t)+A m -KP(t)) lK{t) + A m -KP{t) 



/0 KP+\t) 
where K(t) is a variational parameter satisfying 



(13) 



It is easy to check that, if the predicted closure mean en- 
ergy E*(t) satisfying E*(t) — —A m E P (t) is substituted, 
then r( Gaus;5 )[£; !lt ] = 0. Further insight is obtained by 
considering small perturbations E(t) = E*(t) + SE(t) 
from the predicted mean. By a straightforward calcula- 
tion it follows that 

p(Gcmss) r^j 



8(p - l)A m 
oo fsE(t) + A m ■ pES~ 1 {t)SE{t)) 

, dt - ir^ L+0{SE3) - (15) 



This is the same law of fluctuations as would be realized 
with the Langevin equation 

6E{t) + A rn ■ pET 1 (t)&E(t) = y/2R,(t)r)(t) (16) 

obtained by linearization of the energy-decay equation 
around its solution E* (t) and by addition of a white-noise 
random force rj(t) with strength 



R*(t) 



2(p-l) 



e*(t)£7,(t). 



(17) 



A m -K p (t)+E(t) 



Thus, the smaller fluctuations from the ensemble-mean 
value are predicted to decay according to a linearized law, 
similar to the Onsager regression hypothesis for equilib- 
rium fluctuations. Likewise, the expression Eq.(|l7j) is a 
fluctuation- dissipation relation analogous to that in equi- 
librium. These are testable predictions of the Gaussian 
closure. Note that the coefficient (p — 1) in front of the 
action is > as long as m > —3. In fact, m > —1 is 
required to give a finite energy. Thus, for all permissable 
values of m, the approximate action r( Gauss ) satisfies 
realizability. One should be cautioned again that sat- 
isfaction of realizability is only a consistency check and 
cannot guarantee correctness of predictions. Failure of 
realizability, as observed in the 3-mode model, is more 
practically useful, although in a purely negative way. 

The previous examples and our variational method are 
discussed in greater detail in forthcoming papers Jl5],[l8| . 
Here, we have simply wished to illustrate briefly the use 
of the action principle. Future work will study the suc- 
cess of the new realizability conditions in detecting poor 
closure predictions for more realistic Navier-Stokes flows, 
of greater interest to practical engineering. It should be 
clear that very general PDF ansatz may be employed in 
our method, either by guessing a functional form of the 
PDF or by hypothesizing "surrogate" random variables 
to model the actual flow realizations. Any guess of the 
turbulence statistics — such as the "synthetic turbulence" 
models of ]l9| — may be input to yield predictions for re- 
alistic problems. We therefore expect our method to be 
a flexible framework within which to develop novel tur- 
bulence closures. Insights from simulation, experiment 
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and recent theoretical developments can be readily incor- 
porated. The advantage of the variational formulation 
is that it provides built-in checks of statistical closures 
which may detect a sizable fraction of faulty predictions 
in advance. By doing so cheaply, it can provide great 
savings in turbulence modelling for practical engineering 
purposes. 
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FIGURE CAPTIONS 

Figure 1.) Effective potential for energy in mode 1 in 
quasinormal closure. 

Figure 2.) Effective potential for energy in mode 2 in 
quasinormal closure. (DNS with errorbars). 

Figure 3.) Effective potential for triple moment in 
quasinormal closure. (DNS with errorbars). 
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